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Abstract. The total contribution of strange quarks to the intrinsic spin of the nucleon can be 
determined from a measurement of the strange-quark contribution to the nucleon's elastic axial 
form factor. We have studied the strangeness contribution to the elastic vector and axial form 
factors of the nucleon, using all available elastic electroweak scattering data. Specifically, we 
combine elastic up and up scattering cross section data from the Brookhaven E734 experiment 
with elastic ep and quasi-elastic ed and e- 4 He scattering parity-violating asymmetry data from 
the SAMPLE, HAPPEx, GO and PVA4 experiments. We have not only determined these form 
factors at individual values of momentum-transfer (Q 2 ), as has been done recently, but also 
have fit the Q 2 -dependence of these form factors using simple functional forms. We present the 
results of these fits using existing data, along with some expectations of how our knowledge 
of these form factors can be improved with data from the MicroBooNE experiment planned at 
Fermilab. 



1. The continuing question of the strange quark contribution to the intrinsic spin 
of the nucleon 

The techniques of inclusive and semi-inclusive polarized deep-inelastic scattering employed at 
CERN, SLAC, DESY, and Jefferson Lab have provided a wealth of information about the spin 
structure of the nucleon over the last 25 years. The contributions of the u and d quarks in 
the valence region have now been firmly established. As well, data from collisions of polarized 
protons at RHIC has done much to advance our knowledge of the limitations of the gluon spin 
contribution to the spin of the nucleon. The strange quark contribution to the spin of the 
nucleon, however, is still the subject of investigation, and the indications from deep-inelastic 
scattering are unclear at the moment. 

Consider the determination of As + As from inclusive deep inelastic scattering combined 
with hyperon /3-decay data; a good example would be the HERMES measurement [1] of 
longitudinal spin asymmetries in inclusive positron-proton and positron-deuteron deep-inelastic 
scattering to determine g\ of the proton, deuteron and neutron. This measurement covered the 
kinematic range 0.0041 < x < 0.9, 0.18 GeV 2 < Q 2 < 20 GeV 2 , with the data evolved to Q 2 
= 5 GeV 2 for analysis. Using SU(3) flavor symmetry, these data are combined with triplet and 
octet axial charges (F and D from hyperon /3-decay data) in a NNLO analysis to obtain the 
singlet axial charge and the quark contributions to the proton spin. In doing so, it is necessary 



1 Present address: Department of Physics and Astronomy, Valparaiso University, Valparaiso IN 46383, USA 



to extrapolate the results into the unmeasured regions of x to fill the interval < x < 1. The 
result for As + As, 



As + As = -0.085 ± 0.013(th) ± 0.008(ex) ± 0.009(ev) (HERMES inclusive) 

is inconsistent with to more than 4 standard deviations. 

A different technique is used by the same experiment using data from semi- inclusive deep- 
inelastic scattering, observing asymmetries in the production of charged pions from protons, 
and in production of charged pions and kaons from deuterons [2]. The goal of this analysis is 
to determine the polarized parton distribution functions Aq(x) (and their integrals) over the 
measured x-range only; no extrapolations are performed. As a result, this analysis does not 
rely on SU(3) flavor symmetry for combination with triplet and octet axial charges. However, 
it is necessary to have some understanding the of the fragmentation functions which relate the 
fundamental lepton-quark interaction to the particles observed in the final state. The result for 
As(x) is consistent with zero in the measured x-range; the integral of As(x), 

/■0.30 

/ As(x)dx = +0.028 ± 0.033(stat) ± 0.009(sys) (HERMES semi - inclusive) 

JO.023 

is then of course consistent with zero as well. 

A similar contrasting picture is illustrated by a recent global QCD fit by de Florian, 
Sassot, Stratmann and Vogelsang |3J which is able to bring together the hyperon /3-decay data, 
the inclusive and semi- inclusive deep-inelastic data from CERN, SLAC, DESY, and Jefferson 
Lab, and the pp collision data from RHIC under one roof. Of special interest here is their 
assumption that the strange and anti-strange polarized distributions are equal; this is supported 
by recent COMPASS results [4J. Their functional forms allow for a node in the polarized 
parton distribution functions; this permits a small integral of the strangeness polarized parton 
distribution function if there is cancelation of positive and negative contributions from different 
x-regions. Also, their functional forms allow for SU(2) and SU(3) symmetry violation; however, 
the best fit does not support any significant deviation from these symmetries. They determine 
a "truncated first moment" of the strange quark polarized parton distribution function, 

rl.O 

/ As(x)dx = — 0.006±g'g§f (DSSV truncated) 
io.ooi 

where the uncertainty represents a deviation of the x 2 horn the minimum of the fit by 2%. This 
is clearly consistent with 0. On the other hand, when the full x-range < x < 1 is used, the 
effect of SU(3) symmetry and hyperon b-decay data is seen: 




As{x)dx = -0.057 (DSSV full) 



which implies As + As = —0.114. (The authors of Ref. [3] declined to quote an uncertainty for 
this result because of the uncertainty in the functions required for extrapolation to x = 0.) 

Additional data, of increased precision, from COMPASS at CERN on semi- inclusive deep- 
inelastic scattering [4] only deepens the contrast. Compared to HERMES, this experiment uses 
a rather different polarized beam (muons instead of electrons/positrons), polarized target (solid 
target instead of a gas target), and detector system, but the conclusion reached is similar: 

/■0.3 

/ As(x)dx = -0.01 ± O.Ol(stat) ± 0.01(sys) (COMPASS semi - inclusive) 

J0.004 

As(x) and its integral are consistent with zero in the measured x-rangeJl 
2 See also the contribution by Roland Windmolders. 
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Figure 1. Overview of the types of electroweak elastic scattering data, and their Q 2 values, 
used in the analyses for individual values of Q 2 summarized in Figure [21 as well as in the 
global fit described herein. The data are from the BNL E734 [5J, GO [7], SAMPLE (S], 
HAPPEx Uni QU E2] and PVA4 [TJ El H5] experiments. 

2. Strange quark contribution to the vector and axial form factors of the proton 

It is clearly of interest to examine the strange quark contribution to the nucleon spin in a way 
that is independent of SU(3) symmetry and fragmentation functions. 

The full strange quark contribution to the proton spin, AS, can be directly determined by 
a measurement of the strange contribution to the proton elastic axial form factor, G S A , in low 
energy electroweak elastic scattering. 

AS = As + As = G S A (Q 2 = 0) 

By combining cross sections for vp and vp elastic scattering with parity-violating asymmetries 
observed in eN elastic scattering, the strange quark contributions to the nucleon electromagnetic 
and axial form factors G S E , G S M , and G S A may be determined simultaneously [161 EZ]- The 
reasoning behind this method is briefly reviewed here. 

In elastic scattering of polarized electrons from unpolarized nucleons, one may observe 
a parity- violating helicity-asymmetry in the scattering cross section, A^ v . Combining this 



information with the existing data on the nucleon vector form factors (G E , G E , G P M , G 7 ^), one 
may extract the strange quark contributions G S E and G S M . It is assumed that charge symmetry 
is valid, and that the strange quark distributions in the proton and neutron are identical. This 
idea is the basis of the physics programs of the HAPPEx, GO, SAMPLE, and PVA4 experiments. 
An overview of the data from these experiments is given in Figure [TJ 

But the data on A£ are largely insensitive to the strange quark contribution to the axial 
form factor. To obtain information on G S A , cross sections for vp and Up elastic scattering are 
needed; these cross sections are dominated by the axial form factor at low Q 2 . The only existing 
data for these cross sections are from the BNL E734 experiment [5]. 

A number of analyses have already determined G S E , G S M at individual values of Q 2 , and one 
determines G S A as well. These results are reviewed in Figure [2J Each different analysis has taken 
some subset of data covering a narrow range Q 2 , and uniquely determined the form factors based 
on standard model expressions for the asymmetries and cross sections. The main message from 
these results is that the strange quark contribution to the vector form factors is consistent with 
zero across the full range 0.1 GeV 2 < Q 2 < 1.0 GeV 2 . On the other hand, there is some hint of 
a signal of a negative G S A , but a lack of neutrino-scattering data at low Q 2 prevents a definite 
conclusion about AS = As + As = G S A (Q 2 = 0) at this time. 
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Figure 2. Results of the determination of G S E , G S M , and G S A at individual values of Q 2 , and also 
from our global fit. The separate determinations were done by Liu et al. [18J (green squares at 
0.1 GeV 2 ), Android et al. [7] (blue inverted triangles), Baunack et al. [32] (red squares at 0.23 
GeV 2 ), and Pate et al. [16] (open and closed circles). The preliminary results of the best global 
fit (see text) are shown by the solid line; the 70% confidence level limit curves for the fit are 
shown as the dashed line in the right-hand panel. The vertical scale for G S A in the right-hand 
panel has been adjusted to accommodate the limit curves of the fit. 



3. Global fit of electroweak elastic scattering data 

We have performed a global fit of the available electroweak elastic scattering data, using the 
same general technique described in Ref. Q.6J, but now including all of the data in Figure [1] and 
assuming functional forms for the form factors. A fit of this type can quantity the amount of 
information that can be extracted from the electroweak data in Figure [IJ provide a mechanism 
for including new data that may become available, and can be used as a tool to estimate the 
impact of any future experiments seeking to improve our knowledge of the strangeness form 
factors. 

Our first attempt to fit the form factors G E , G S M , and G S A used this simple set of functional 
forms: 

G E= (l + Q2 T /A|)2 Gll= (l + QV A 2 f)2 G A= {1 + q2 /A 2 a) 2 (initial function set) 

where r = Q 2 /4M^, p s = (dG s E / 'dr)| T= o is the strangeness radius, /i s is the strangeness magnetic 
moment, and the As are shape parameters that define how rapidly these functions depend on 
Q 2 . These functions satisfy a number of boundary conditions: they all depend only on even 
powers of Q; they tend towards as Q 2 — > oo, and they come to the correct physical value at 
Q 2 = 0. We found, however, that this set of functions is not quite suited to the existing data. 
First of all, a good fit to the neutrino data was not possible unless very large negative (and 
clearly unphysical) values of AS were permitted; this problem was solved by introducing a term 
linear in Q 2 in the numerator of the G S A function. Secondly, the shape parameters A^ and Am 
are not determined; the best fit values for them are very large, with very large errors. This is 
not surprising; looking at Figure [2] it seems that both G S E and G S M are essentially featureless, 
since they are consistent with zero. If A E and Am are very large, then the Q 2 -dependence in 
the denominators of G S E and G S M may be neglected and we may use a simple linear function for 
G S E and a constant function for G S M . Then our final set of functions were these. 

G% = PsT G S M = p s G A = ^j-^ q 2 /A^)2 ( final functi ° n Set ) 

The best values for the five parameters (preliminary) are: 

p s = 0.13 ±0.21 fi s = 0.035 ± 0.053 AS = -0.27 ± 0.41 A j4 = 1.3±1.9 S A = 0.32 ± 0.48 

The best fit is shown as the solid line in Figure [2J the 70% confidence level uncertainty limits 
are shown by the dashed line. 

It should be stressed that the fit has not been made to the values of the form factors displayed 
in Figure [2j Instead, the fit has been made to the measured values of the cross sections and 
asymmetries from the experiments mentioned in Figure [TJ In Figure [21 we show that this global 
fit is consistent with the results of previous limited analyses of subsets of these data. 

The preliminary results of the fit can be simply described. The strangeness radius and 
magnetic moment are consistent with zero, and the uncertainties in these parameters are 
consistent with the uncertainties in the separate determinations of G S E and G S M on display 
in Figure [2j On the other hand, AS is also consistent with but the uncertainty is very large 
because there are no vp or vp elastic data at sufficiently low Q 2 to constrain it. As a result the 
uncertainties in the global fit to G S A are very much larger than the uncertainties in the separate 
determinations of G S A in Figure [2j 

It is clear that the single greatest need is for a much improved set of data on neutral current 
vp and vp elastic scattering down to the lowest Q 2 possible. In the last few years it has become 
clear that a better understanding of the charged-current cross section is also needed, as described 
in the next section. 



4. The problem of the charged-current axial form factor 

The results of this fit are still preliminary because we do not have a full systematic uncertainty 
analysis. Such an analysis would consider variations in all the fixed input parameters, such as 
the full nucleon form factors, within their known limits as determined by existing experimental 
data or by theoretical considerations. 

One nucleon form factor that must be known for this analysis to proceed is the charged-current 
portion of the axial form factor. The full nucleon axial form factor contains contributions from 
up, down, and strange quarks, 

G A = - (-G U A + G d A + G s A ^j 

and this form factor occurs in the current associated with Z-exchange processes like ep, up, and 
vp elastic scattering. The "charged-current" portion of that form factor, 

A — ~ <-M 

occurs in W-exchange processes, like Vu_ + n — > [i~ + p. The cross section for i/„ + n — > \i~ + p 
has been studied using deuterium and nuclear targets for many years, and until a few years 
ago the form factor G A extracted from those cross sections was considered "known." In this 
analysis we used the parametrization provided by Budd, Bodek and Arrington [W \ I2UI |2"T] . In 
the analysis technique used here and in Ref. [16], it is the full axial form factor G A which is 
determined from the neutral current ep, up, and up data; the charged-current part G A must 
be subtacted away to obtain G S A , any uncertainty in G A translates directly into an uncertainty 
inG^. 

In recent years, the charged-current portion of the axial form factor has become a lot less 
"known" than it used to be, due to high-statistics neutrino-nucleus scattering experiments such 
as MiniBooNE [22j [23] and K2K [2l] that raise questions about our knowledge of its (In- 
dependence. Due to the imperfect knowledge of the incoming neutrino momentum and inability 
to reconstruct the vertex kinematics, a more complete and sophisticated way of analysing these 
data must be developed which takes into account the electroweak nuclear response. This is 
underway - see for example Benhar, Coletti, and Meloni [25], Ankowski, Behnar, and Farina [26], 
and references therein. 

In order to improve the situation, it is necessary to have both the charged-current and neutral- 
current elastic channels measured at low Q 2 in the same experiment and treated in a consistent 
way to extract the axial form factor (s). The MicroBooNE experiment at Fermilab may provide 
these data. 

MicroBooNE is an approved experiment at Fermilab to build a large liquid Argon Time 
Projection Chamber (LArTPC) to be exposed to the Booster neutrino beam and the NuMI beam 
at Fermilab. The experiment will address the low energy excess observed by the MiniBooNE 
experiment, measure low energy neutrino cross sections, and serve as the necessary next step in 
a phased program towards massive Liquid Argon TPC detectors. 

Liquid argon is very well-suited to observe low energy protons from low-Q 2 neutral-current 
and charged-current scattering; this experiment may be an ideal place to move towards a 
successful determination of As. 

5. Summary 

The strangeness contribution to the nucleon spin may be determined from a combined analysis 
of low-energy electron-nucleon and neutrino-nucleon elastic scattering data. 

However, existing neutral-current neutrino data (from BNL E734) lack insufficient precision 
and Q 2 -range to make possible a definitive determination. In addition, recent neutrino 



experiments have called into question our understanding of the charged-current axial form factor 
extracted from data on nuclear targets. 

New experiments (e.g. MicroBooNE) can provide the datasets needed for a consistent 
treatment of the electroweak nuclear response and extraction of the strange axial form factor. 

Acknowledgments 

This work was funded by the US Department of Energy, Office of Science. 
References 

[1] Airapetian A et al. (HERMES) 2007 Phys. Rev. D75 012007 (Preprint hep-ex/0609039 ) 

[2] Airapetian A et al. (HERMES) 2005 Phys. Rev. D71 012003 (Preprint hep-ex/0407032 ) 

[3] de Florian D, Sassot R, Stratmann M and Vogelsang W 2009 Phys. Rev. D80 034030 (Prevrint [0904. 382 ip 

[4] Alekseev M G et al. (COMPASS) 2010 Phys. Lett. B693 227-235 (Prevrint Il007.406ip 

[5] Ahrens L A et al. 1987 Phys. Rev. D35 785 

[6] Armstrong D S et al. (GO) 2005 Phys. Rev. Lett. 95 092001 ( Preprint |nucl- ex /0506021 ) 
[7] Android D et al. (GO) 2010 Phys. Rev. Lett. 104 012001 (Prevrint 10909 .5107|l 

[8] Beise E J, Pitt M L and Spayde D T 2005 Prog. Part. Nucl. Phys. 54 289-350 ( Preprint |nucl- ex/04120541 
[9] Aniol K A et al. (HAPPEx) 2004 Phys. Rev. C69 065501 ( Preprint |nucl-ex/0402004l 
[10] Aniol K A et al. (HAPPEx) 2006 Phys. Lett. B635 275-279 (Preprint [n"ucl-ex /0506011 ) 
[11] Aniol K A et al. (HAPPEx) 2006 Phys. Rev. Lett. 96 022003 ( Preprint |n"ucl-ex/0506010| ) 
[12] Acha A et al. (HAPPEx) 2007 Phys. Rev. Lett. 98 032301 ( Preprint |n"ucl-ex/0609002T ) 
[13] Maas F E et al. (A4) 2004 Phys. Rev. Lett. 93 022002 ( Preprint [5ucl-ex/0401019l 
[14] Maas F E et al. (A4) 2005 Phys. Rev. Lett. 94 152001 ( Preprint [nucl-ex /0412030 ) 
[15] Baunack S et al. (A4) 2009 Phys. Rev. Lett. 102 151803 (Prevrint l0903~2 733) 
[16] Pate S F, McKee D W and Papavassiliou V 2008 Phys. Rev. C78 015207 (Preprint 0805 . 2889) 
[17] Pate S F 2004 Phys. Rev. Lett. 92 082002 (Preprint |hep- ex/0310052) 

[18] Liu J, McKeown R D and Ramsey-Musolf M J 2007 Phys. Rev. C76 025202 (Prevrint [0706 ■ 0226P 

[19] Budd H, Bodek A and Arrington J 2003 (Preprint [hep- ex/ 0308005 ) 

[20] Bodek A, Budd H and Arrington J 2004 AIP Conf. Proc. 698 148-152 

[21] Budd H, Bodek A and Arrington J 2005 Nucl. Phys. Proc. Suppl. 139 90-95 

[22] Aguilar-Arevalo A A et al. (MiniBooNE) 2008 Phys. Rev. Lett. 100 032301 (Prevrint 10706.09260 
[23] Aguilar-Arevalo A A et al. (MiniBooNE) 2010 Phys. Rev. D81 092005 (Preprint 1002 . 2680) 
[24] Gran R et al. (K2K) 2006 Phys. Rev. D74 052002 (Preprint [hep- ex/ 0603034 ) 
[25] Benhar O, Coletti P and Meloni D 2010 Phys. Rev. Lett. 105 132301 (Preprint [T606. 4783 P 
[26] Ankowski A M, Benhar O and Farina N 2010 Phys. Rev. D82 013002 (Prevrint [TOOT . 0481) 



